Sommerfeld's image method in the calculation of van der Waals forces 
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We show how the image method can be used together with a recent method developed by C. 
Eberlein and R. Zietal to obtain the dispersive van der Waals interaction between an atom and a 
perfectly conducting surface of arbitrary shape. We discuss in detail the case of an atom and a semi- 
infinite conducting plane. In order to employ the above procedure to this problem it is necessary to 
use the ingenious image method introduced by Sommerfeld more than one century ago, which is a 
generalization of the standard procedure. Finally, we briefly discuss other interesting situations that 
can also be treated by the joint use of Sommerfeld's image technique and Eberlein-Zietal method. 

PACS numbers: 



I. INTRODUCTION 

Recent technological advances are responsible for the 
increasing number of experiments in the micro and nano 
scales. In such scales, dispersive forces become very im- 
portant and, consequently, much effort has been made 
in better understanding the role played by these forces. 
Many papers have been dedicated recently to the problem 
of atoms interacting with different surfaces, see 
and references therein, to mention just a few. 

In this paper we show how to employ a recent method 
developed by C. Eberlein and R. Zietal Q in connec- 
tion with the image technique largely used in electrostatic 
problems in order to obtain analytically expressions for 
the non-retarded dispersive interaction energy between 
an atom and a perfectly conducting surface. Particu- 
larly, we explore this connection to obtain the dispersive 
van der Waals interaction between an atom and a semi- 
infinite conducting plane. Although one does not expect 
that the image method can be applied to this problem, 
an ingenious extension of the image method, provided by 
A. Sommerfeld in the late nineteenth century, enables us 
to treat this and other interesting non-trivial systems by 
such a method. 

This paper is organized as follows: in section 2 we 
briefly review Eberlein-Zietal procedure and show how 
the image method can be used with it to yield the non- 
retarded dispersive interaction between an atom and a 
perfectly conducting surface. In section 3 we show how 
to apply Sommerfeld's image method to obtain the elec- 
trostatic potential in the non-trivial system constituted 
by a point charge in the presence of a semi-infinite con- 
ducting plane. Since Sommerfeld's work is largely for- 
gotten, we dedicate ourselves, in this section, to present 
a pedagogical exposition of Sommerfeld's main ideas. In 
section 4 we show how the methods introduced in sec- 
tions 2 and 3 can be used together to yield the quantum 
(non-retarded) dispersive interaction between an atom 
and a semi-infinite conducting plane, a result previously 
computed by other method by Eberlein and Zietal Q • In 
section 5 we discuss other geometries where Sommerfeld's 
image technique is also useful, giving special attention to 
the atom-disk system and the system constituted by an 



atom and an infinite plane with a circular aperture. Sec- 
tion 6 is left for the final remarks. 



II. IMAGE METHOD AND THE VAN DER 
WAALS FORCE 

The use of the image method has rendered fruitful 
results in electrostatics, hydrodynamics, heat equation, 
wave equation, etc. @~[l2|. A recent method developed 
by C. Eberlein and R. Zietal Q brings the possibility of 
a very simple and natural way of employing the image 
method for the calculation of the van der Waals disper- 
sive interaction of an atom and a perfectly conducting 
surface of arbitrary shape. In this section we briefly 
state this method and illustrate it with a very simple 
and well-known system, namely, an atom near an infinite 
conducting plane. 

Let a neutral atom be at position ro in the presence of 
an arbitrary perfectly conducting surface S. In the non- 
retarded regime, i.e., for distances much shorter than the 
dominant atomic transition wavelength, the dispersive 
interaction energy can be written in the form Q 
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where di is the i component of the atomic dipole opera- 
tor d and (•••) means quantum expectation value taken 
with the atom in its ground state. The above expression 
is valid for any orthonormal coordinate system and for 
atoms with no permanent dipole moment. Both condi- 
tions can be easily abandoned with only minor changes 
in the formula. For grounded conductors, the function 
Gh appearing in equation ([!]) satisfies the equation 



V 2 G H (r,r') =0, 
submitted to the boundary condition 



f 



47r|r — r'| 



= 0. 



(2) 
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For non-grounded conductors, the boundary condition 
must be modified. Previous equations are very similar 
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to those satisfied by the potential Vi (r) generated by the 
image charges in the electrostatic problem of a charge q 
at position r' in the presence of a perfectly conducting 
grounded surface S, if the geometry accepts an image 
treatment. Vi satisfies Laplace equation and a boundary 
condition that differs from ([3]) only by a constant factor. 
We are thus allowed to make the identification 



G H (r,r') 



'1 



(4) 



Therefore, to hnd the dispersive van der Waals interac- 
tion between an atom and an arbitrary perfectly con- 
ducting surface all one ultimately needs is to solve an 
electrostatic problem. Afterwards, by employing formula 
tH), one can obtain the desired quantum dispersive van 
der Waals interaction in a straightforward way. This is 
not surprising since in the short distance regime (non- 
retarded regime) the retardation of the fields can be ne- 
glected, allowing us to remain in the static regime. For 
distances of the order of the wavelength of atomic tran- 
sitions, where retardation effects become important, it 
is also possible to change from a quantum problem to 
a classical one, but it demands to solve an electrody- 
namic problem, dealing with Hclmholtz equation instead 
of Laplace e quat ion. This is the essence of the scattering 
method, see [13[ and references therein. For pedagogical 
reasons, we close this section by illustrating the proce- 
dure just described in a very elementary situation. Let 
the surface S be an infinite conducting plane located at 
y = 0. The image method corresponding to this geome- 
try is trivial and leads us, with the aid of (@|, to 



G H (r,r') 
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(5) 

Substituting this solution into cq. ([I]) we obtain the well- 
known result [lJI 
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This example shows the simplicity of the image method in 
dealing with dispersive interactions. The pitfall seems to 
be the small number of geometries compatible with this 
kind of treatment. Nevertheless, that number was fortu- 
nately increased with the ingenious extension developed 
by Arnold Sommcrfcld in the end of the 19th century, as 
we shall see along this work. 



III. SOMMERFELD'S METHOD FOR A 
CHARGE AND A SEMI-INFINITE PLANE 

To begin with, consider a charge q in the presence of a 
semi-infinite perfectly conducting plane, taken as y = 0, 
x > 0. As discussed in the last section, to find the elec- 
trostatic potential at every point of the space (where it 
is defined) for this problem is all we need to obtain the 



van der Waals interaction between an atom and a semi- 
infinite conducting plane. The appropriate system of co- 
ordinates for this problem is the cylindrical set (r, tp, z), 
since with this choice the equation of the semi-plane is 
only tp = 0, 2-7T. Naively, one may assume that image 
method is unsuited for this geometry once all points of 
space seem to be already used, leaving no room for im- 
age charges. Sommerfeld's procedure [lj|'[l6[ is based 
on creating somehow space to allocate the image charges. 
In his construction, to come back to the point of depar- 
ture a rotation of 4-zr is required. In other words, the 
points (r, tp, z) and (r, tp + 2n, z) are no longer the same. 
In this way, Sommcrfcld revived the riemannian sheets 
of complex analysis, by constructing a two-fold space - 
points with < tp < 2tt belong to the ordinary space 
while those with 2ir < tp < Att constitute the auxiliary 
space, where we will be able to locate the necessary im- 
age charges. Note that with this construction we create 
a discontinuity at the conducting surface. Indeed, when 
we approach the semi-plane from above in the ordinary 
space we have tp — > 0, while when we approach from 
below, we have tp — > 2tt. Let us consider the charge q 
at position (r',tp',z'). The potential generated by this 
charge in the double space cannot be the familiar 1/R, 
with R = |r — r'|. This is not surprising since the poten- 
tial created by a point charge depends on the geometry 
of the space - recall that in 2 dimensions the coulomb po- 
tential is not any more proportional to 1 over the distance 
from the point of space and the position of the charge, 
but it is given by a logarithm of such a distance over a 
reference distance. The potential 1/R presents a symme- 
try by changing tp by tp+2n, so its laplacian furnishes two 
Dirac delta functions, one with singularity at (r',tp',z') 
and another with singularity at (/, tp' + 2ir, z'). We are 
thus led to recognize in 1/R the superposition of the po- 
tential of two distinct charges, one in the ordinary space 
and the other in the correspondent place of the auxiliary 
space. In order to identify each term, we employ Cauchy 
theorem, which states that 



R(w) 
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2ttI 



dw . 
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(7) 



where C is any closed contour in the complex plane, pro- 
vided that it encloses the pole w' = w and R~ 1 (w') is 
analytical in the interior of and along C. Wc must chose 
a convenient variable w. Since tp is the variable which 
operates the passage from the ordinary to the auxiliary 
space, the smart choice is 

i V /2 



W 



(8) 



which has a period of in, as the double space. Accord- 
ingly, we set w' = e la l 2 which leaves equation in the 
form 
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1 



1 



,ia/2 



R{(p) 4tt Jc R{a) e la / 2 - e^l 



-da . 



(9) 



Now, all we have to do is to choose the contour C. The 
distance between points (r,tp,z) and (r' ,tp' , z') is given 
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FIG. 1: Choice of the contour, 
the dashed lines. 



The cuts are represented by 



by R 2 = r 2 + r' 2 - 2rr' cos(p - tp') + (z - z') 2 . R(a) is 
obtained from R just changing p by a. The presence of 
the square root in R^ 1 (a) generates a branch cut. The 
branch points are determined by the equation R(a) = 0, 
whose solutions are a = ip' + 2m7r ±17, mGZ, with the 
definition 



cos(i7) = 



+ r' 2 + {z- z' f 
2rr' 



(10) 



Our contour must avoid the cut generated on these 
branch points. One possibility is sketched in the Fig- 
ure [U Now, we shall perform the integration ([9]) in this 
contour. Due to the periodicity of the integrand, the in- 
tegrations at vertical lines Re a = and Re a = An cancel 
out. The contributions of the horizontal paths vanish in 
the limit where the height of the circuit in Figure [1] is 
taken to infinite. This is evident once i?" 1 — > in this 
limit, leaving the integrand of <j9j> null there. We are left 
with the integrals around the cuts. Let us call A the 
path around the branch points ip' ± 17 and A\ the one 
around the branch points <p' + 2n + ±17. Therefore, the 
potential 1/R may be written as 
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4tt 



R,, 
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da+ — I ^ —da. 
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(11) 

We achieved a separation of the potential 1/R in two 
terms. Sommerfeld showed that the first (i) is uniquely 
defined, finite and continuous at all points of the dou- 
ble space, except at (r' ,cp',z'). This means, particularly, 
that it is finite at (r', p' + 2tt, z'); (ii) has a null laplacian 
at all points except at (r', ip' , z') and at the surface of the 
conductor; (Hi) vanishes at infinity and (iv) is bivalent in 
the ordinary space, with a separated branch to each copy 
of the double space. Therefore, including the necessary 
constants, we realize that the potential of a single charge 



at position (r', ip' z') in the double space is identified as 



V(r,(p,z) 



<1 



16tt 2 £o J 1 - e^-")/2 
^0 



da . 



(12) 



Analogously, the term containing the integration along 
the path A\ can be identified as the potential in the dou- 
ble space created by a charge at position (/, ip' + 2ir, z'). 
To evaluate the integral in the previous equation, we 
shall first observe that the contour A$ has two paths, one 
around each branch point, and each path contains three 
parts, two verticals and a semi-circle around the branch 
point. In the limit where the radius of this semi-circle 
goes to zero, the integral along the semi-circle vanishes. 
The cut prevents the contributions of the two vertical 
lines to cancel each other. It can be shown from complex 
analysis that they give the same result. In order to calcu- 
late it, we parametrize the vertical line around the branch 
point a — ip' + i"f by a = ip' + i/3, and the one around the 
branch point below, a = ip' — i'y, by a = ip' — i/3. This 
way, both integrals run with varying from 7 to 00 , and 
a factor 2 must be included to account for both of them. 
Hence, equation (fl~2"|) becomes 



V(J>,6,<P) = 



8TT 2 Eo\ / rF 

(cosh P - cosh7)" 1 / 2 sinh(^/2) 



cosh(/3/2) 



v-v' 

2 



dp. (13) 



This integral can be recast with a convenient change of 
variables. Defining 



£ := cosh(/3/2) , a := cosh(7/2) , r := cos 



ip-ip 



we rewrite integral (fT3j) as 



(14) 



V(r,ip,z) 



<l 



8TT 2 e V2r 



r/£ 



(£2_ a2) l/ 2( £_ r) 



q 



2TT 2 e R 



tan 



(a-r) 



1/2 



(15) 



Now that we know the potential of a point charge in the 
double space, we may proceed with the image technique 
to find the electrostatic potential of a charge q in the 
presence of a semi-infinite plane. It is easily seen that if 
we put an image charge —q at (r',47r — p',z') we have 
the boundary conditions satisfied. Indeed, the potential 
of the system constituted by these two charges is 



V sp (r,(p,z) 



2ir 2 e R 

q 

2n 2 £ R t 



tan 



1/2 



tan 



{a — 1 

(<r + Tj) 
(cr - n) 



1/2 



(16) 



where Rj and Tj are obtained from R and r just chang- 
ing ip' by 47r — ip' and a is the same for both charges 
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since it doesn't depend on ip' . At the conducting surface, 
characterized by ip = 0, 2-7T, we have R = Ri and t = t%. 
Therefore, V sp vanishes at the conducting semi-infinite 
plane. Observe that for <p = n, we have V sp 7^ 0, as 
expected. Note, also, that the position occupied by the 
image charge is the correspondent point in the auxiliary 
space of the mirror image of the charge. Henceforth, we 
recognize the potential (fTo| as the potential of the con- 
figuration formed by a charge and a semi-infinite plane. 



IV. AN ATOM AND A SEMI-INFINITE PLANE 

To employ relation ([4]) we must first determine the 
potential Vi created by the image charge in the chargc- 
semi-infinite plane configuration. This is easily done once 
we know the complete potential, given by (|16p . 

<1 



Vi(r,ip,z) = V sp (r,(p,z) - ^ 
Substituting this into equation ([4]) we obtain 



(17) 



Gh(t,t') = - 



1 



2tt 2 R 
1 



tan 



tan 



1/2 



nl/2 



(18) 



Applying (JTJ) , we obtain the desired van der Waals inter 
action between an atom at ro and a semi-infinite plane, 

-1 



V. OTHER NON-TRIVIAL APPLICATIONS 

As already mentioned, analytical expressions of disper- 
sive interactions for different geometries are important in 
a variety of situations. Hence, the natural question one 
may pose is whether there are other geometries which 
admit a treatment by the image method. The purpose 
of this section is to answer this question positively. In 
the solution for the semi-plane, we employed a two-fold 
space. Some times, we need to deal with more spaces. 
As an example, think of an atom in the presence of a 
wedge. If the angle of the wedge is n/m, with m being 
a positive integer, we may employ the standard image 
technique to solve the electrostatic problem and then ap- 
ply Eberlein-Zietal formula. On the other hand, let the 
angle be nir/m, with n, m positive integers. The stan- 
dard image procedure yields images on the region inside 
the wedge, leading to the failure of the method. In spite 
of this, we may approach the problem through Sommer- 
feld's generalization. Consider, without loss of generality, 
an irreducible fraction n/m. To apply the formalism of 
the last section we will have to deal with n spaces. This 
makes the calculations a little harder, but the ideas are 
the same. To give a sample, the fundamental modifica- 
tion is the variable we choose in Cauchy's theorem (|7J). 
Since the n-fold space has a period 27m, instead of ([8]). 
we must choose 



W 



Jtp/n 



(22) 



47T£ 



[A(dj 



B(d%) 



+ C(d 2 z )} , (19) As a consequence, we must substitute equation ([9]) by 



with 
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COS (fo 
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<A)) 



= i<M- 1 («)- 



R(tp) 27m 



v/v; 



da . 



(23) 
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167rrQ sin tpo 

cos ipo 



167rrg sin ipo 
|tt - (f \(l + cos 2 ip ) 



C 



1 



247rr3 



8nr% sin 
cos ifo 
167rrn sin 2 



Vo 



167T7q sin 

■ ~ Vol 



167rrn sin 



(20) 

Note that, by symmetry arguments, the interaction en- 
ergy is independent of Zq. This problem was analysed 
by C. Eberlein and R.Zietal Q who solved it by employ- 
ing properties of Bessel functions. It is possible to show 
that the Gh obtained there for this geometry is equiv- 
alent to our equation ([T8|) . It is instructive to analyse 
our results in the limit where the atom is very close to 
the conducting surface and far from the edge. To do 
so, we must take the limits r — >• 00, <p a — >• 0, keeping 
the product 7/0 = r sm(p constant. Doing this, we get 

A -> l/(16yo) i B -> 1 /i 8 Vo) and c ~* 1 /i 16 Vo), which 
allows us to recast the interaction energy in the form 

+ 2(%) + {<%)], (21^ 



This potential is the superposition of the potential of n 
point charges, located at (r', p'+i2n, z'), (i = 0, .., n— 1). 
The decomposition of the potential l/R into n contri- 
butions is done in a similar fashion. We must enlarge 
the circuit of Figure [1] by changing the vertical line at 
Re a = Att by the vertical line at Re a = 2irn, so that af- 
ter integration the vertical contributions cancel out due 
to the periodicity of the integrand. This enlargement 
brings more cuts to the integral, leading to the desired 
extra terms in (jll[) . which becomes 



1 

R 



" ± 1 

Y — 

t— 1 2nn 

3=0 A 



R„ 



l — e i(v-ai)/2 



da . 



(24) 
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As expected, last result coincides with that obtained for 
an atom in the presence of an infinite plane, see equation 
©• 



From the operational point of view, the only change from 
the two-fold space to the n-fold one is the integration we 
must perform. For more details, see the original memoir 

In the last paragraph we saw other possibilities of em- 
ploying the image technique increasing the number of 
folds of the space. We may find also examples that re- 
quire only the two-fold space treated before. To this end, 
we generally have to use different system of coordinates. 



5 



It is of foremost importance in the image method to de- 
scribe the problem by the system of coordinates that suits 
it. A relevant example is that of an atom interacting with 
an infinite plane with a circular aperture, which leads to 
repulsion as shown by Levin et al.0. In a recent pa- 
per [I?]], we treated this problem by Sommerfeld's image 
method and obtained an analytical result for the corre- 
sponding interaction energy of this system. Our results 
agree with those previously obtained by Eberlein and Zi- 
etal H. We close this section by summing up our main 
results for a system constituted by an atom and an in- 
finite conducting plane with a circular aperture and for 
the atom-disk system. 

The interaction energy between an atom and an infinite 
plane with a circular aperture is found in a similar way 
as we did in the semi-plane system. For an atom on the 
symmetry axis of the aperture, say the z-axis, and with 
major pollarizability in the z direction, our expression 
becomes very simple and it is given by 



E p h = 



(4 



1 H — sin 

7T 



4az(3a 4 + 8a 2 z 2 - 3z 4 ) 
37r(a 2 + z 2 ) 3 



(25) 



where a is the radius of the aperture. For small z 
we get repulsion. We found an equilibrium point at 
z e q ~ 0.74235a. For z > z eq the interaction is attractive 
so that the equilibrium is stable in the z direction. How- 
ever, it can be shown that for lateral displacements the 
equilibrium is unstable. Atomic anisotropy is essential to 
get repulsion. With this method we may also treat the 
complementary geometry, namely, the atom-disk system. 
The interaction energy when the atom is in the symme- 
try axis of the disk and with dominant pollarizability in 
this direction is [13 



E, 



d isk 



(4 



1 sin 

7T 



+ 



64e 7rz 3 

4az(3a 4 + 4a 2 z 2 + 9z 4 ) 
37r(a 2 + z 2 ) 3 



z 2 + a z 



(26) 



where a is now the radius of the disk. Expressions 
(|25"j) and (|26|) allow us to make a quantitative study 
of the finite-size effect for those geometries. Other in- 
teresting feature is the possibility of an exact calcula- 
tion of the non-additivity effects. Indeed, evaluating 
Fdisk = -d z E dlsk and summing it with F ph = -d z E ph 
we obtain 



F, 



disk 



F, 



ph 



F n 



(d 2 z )a (z 2 - a 2 )z 



e ir 2 (z 2 + a 2 ) 4 



(27) 



where Fq is the atom-infinite plane force. We see that 
non-addivity vanishes for a — > oo and for a — > 0, as it 
should. It is remarkable the existence of a distance for 
which the superposition holds. In eq. ([2T|) this distance is 
found to be z = a. 

VI. FINAL REMARKS 



The main purpose of this paper was to establish the im- 
age method as a powerful tool to calculate non-retarded 
dispersive forces between an atom and a conducting sur- 
face. The recent procedure developed by Eberlein and 
Zietal is remarkably and naturally linked to the elec- 
trostatic image method. As an example, we discussed in 
detail the non-trivial system constituted by an atom and 
a semi-infinite conducting plane. This problem reveals 
the importance of choosing an adequate system of coor- 
dinates, since it is the angular variable ip which sustains 
and operates the two-fold construction of the Sommer- 
feld's image method. We briefly analysed in a similar 
way two other non-trivial systems, namely, an atom and 
an infinite conducting plate with a circular aperture and 
the atom-disk system. Sommerfeld's image method is not 
systematic in the sense that for each different geometry 
one has to work out the proper system of coordinates 
and develop the multi-fold space. On one hand, this 
may be considered a pitfall but, on the other hand, it 
leaves room for the creativity and imagination of the re- 
searcher when challenged to solve non-trivial problems. 
Once the exact analytical dispersive interaction energy 
is obtained, a variety of questions may be answered. In 
this connection, we would like to mention particularly 
the possibility of studying finite size effects and, in some 
cases, non-additivity effects involving complcmctary sur- 
faces. In the latter case, the complementary geometries 
of an atom-disk system and the system constituted by an 
atom and a plane with a circular aperture was studied in 
[[13] i where a surprising result was found, namely: that 
there exists a distance between the atom and the center 
of the aperture for which the non-additivity effects van- 
ish. The authors believe this result may occur in other 
systems involving complementary surfaces. 
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